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Benchmark helium dimer and trimer calculations with a public 

few-body code. 

Vladimir Roudnev and Michael Cavagnero 
Abstract 

We present detailed calculations of bound and scattering states of dimers and trimers of He to 
produce highly accurate data and to test a non-relativistic three-body code currently in develop- 
ment for public distribution. For these systems, uncertainties and inaccuracies in the fundamental 
constants frequently used in published works can substantially exceed numerical errors. Our bench- 



i-£^ ■ mark calculations include specific estimates of the numerical accuracy of the calculations, and also 

Ph; 

Oj. explore sensitivity to fundamental constants and their uncertainties. The use of an inexact cou- 

pling constant in the previous calculations leads to 0.08% error for the ground state energy, 0.3% 
error for the excited state energy and up to 0.15% error for the atom-dimer scattering length in 
the system of three ^He atoms. The corresponding errors for the unsymmetric ^He2'^He system are 
0.3% for the bound state energy and 0.03% for the atom-dimer scattering length. 



I. INTRODUCTION 

Experimental investigations of He dimers and trimers [l|-l3| in the 1990's generated con- 
siderable theoretical interest in these weakly-bound few-body systems. Small He clusters 
are now considered one of the most extensively studied molecules - theoretically and com- 



putationally 



- |31| (see also an extensive review [31)])- Many semi-empirical and ab initio 



potential models have been proposed for the interaction of He atoms |4J-llO|. These potential 
models have been^used as input in a number of two- three- and even four- body quantum 
calculations 
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26| . The numerical results of these investigations are, for the most part, in 
agreement, with small remaining discrepancies due simply to limits of numerical accuracy. 
Nevertheless, we wish to revisit, reproduce and check these previous results. While small 
clusters of He are structurally quite simple, they are also delicate, with very weak binding, 
and so are challenging to model numerically. Given the extensive set of known potential 
models and the number of results published, it seems natural to use bound and scattering 
states of two and three He atoms to benchmark the accuracy of existing and future few-body 
quantum-chemical software. 

The near-threshold state of the He dimer is extremely sensitive to the details of the 
interaction potential. Visually indistinguishable potentials can lead to qualitatively different 
expectations for the properties of the He dimer. The three-body states are likewise expected 
to be very sensitive to the details of the interaction, since their binding is similar to that of 
Efimov states. One immediate result of our investigations is that, even for a fixed potential, 
He cluster observables can vary substantially depending on the accuracy of the fundamental 
constants used in the calculations. Unfortunately, the fundamental constants and conversion 
factors published by the developers of potential models and early He cluster calculations - 
and, therefore, widely used in subsequent publications - have intrinsic systematic errors that 
prevent them from being used for benchmark purposes. 

In this paper, we consider three questions: To what extent is the accuracy of near- 
threshold calculations limited by knowledge of the fundamental constants? To what extent 
can the known discrepancies in the published results be explained by possible use of inexact 
coupling constants? Can our three-body calculations reach the "natural" accuracy limit set 
by the uncertainties of the known fundamental constants? Answers to these questions will 
allow us to benchmark a quantum few-body code being prepared for public release, with 



intended applications to a wider variety of atomic and molecular few-body systems. 

The detailed technical description of the three-body code we use for the calculations will 
be made elsewhere; gene rally, the techniques used in these calculations are similar to the ones 



employed in 18 
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25| . In the current stage of development the code is not yet completely 
ready for an extensive public distribution. Interested readers, however, are encouraged to 
contact the authors for obtaining the current development version. 

We assume throughout this paper that the system of He atoms can be described by 
a single-channel pairwise interaction. We, therefore ignore the contribution of three-body 
forces. This contribution can, substantially exceed the numerical accuracy of the calculations 
presented. Our goal, however, is to provide an exemplary data set based on known potential 
models rather than to perform the most realistic calculation of the physical system. We 
also ignore all the off-diagonal non-adiabatic corrections not included in the single-channel 
potential model. 

In the following sections we shall give a short description of the problem, analyze the 
uncertainties that enter the benchmarking problem, give a review of published results and 
provide our results, which account for all known physical uncertainties. 

II. PHYSICAL UNCERTAINTIES IN TWO-BODY OBSERVABLES 

Consider the Schrodinger equation for a two-body system 

{-—V^ + cV{r)-E)m = Q . 
2m 

Here, c is a conversion factor between the energy units of the potential (provided by the 
various developers of two-body potentials) and our chosen "natural" units, as specified below. 
This conversion factor, specified with varying degree of precision by different authors, is 
usually a simple function of fundamental constants. Rescaling the energy we reduce the 
equation to 

{-V^ + v{r) - z)"^ = , 

where v(r) = c^V(r) and z = ^E. To analyze the sensitivity of the state \E' to the 
inaccuracy of the fundamental constants and other conversion factors, we introduce a (di- 
mensionless) coupling constant A and treat the eigenvalue z as a function of the coupling 



constant 

i-V"^ + Xv{r) - z{X))^ = . 

Knowing the derivative of the energy eigenvalue ■£^z{\)\x=i one can estimate the sensitiv- 
ity of the dimer bound state, or scattering phase shifts, to the degree of precision of the 
fundamental constants of interest. 

Many model potentials of rare gas atom interactions were published in units of tem- 
perature. The conversion factors between atomic units of energy (typically employed in 
calculations) and temperature units are, however, not always provided. In some cases, when 
the corresponding conversion factors are published, the values suggested by the authors differ 
significantly from the values recommended by NIST. For instance, the recommended value 
for the TTY potential is 3.1669x10"^ a.u. K~^, which differs from the recommended value 
of 3.1668154x10^® a.u. K~^ in the 5th significant figure. Even such small differences can 
appreciably affect the accuracy of derived results owing to the weakly-bound nature of few- 
body clusters, and can make benchmarking and comparison of alternative computational 
strategies difficult. 



Many three-body calculations of He clusters 12|, [iSl, |26|, |3l|] have been performed, as- 



suming A for distance and the conversion coefficient -^^^ ~ 12.12 A^K. Keeping only 4 
significant figures may lead, however, to significant inaccuracy in 2-body observables. Let 
us estimate the accuracy of this coupling constant based on the best available data from the 
NIST database. Based on the following recommended values 

2/i = m4He = (4.0026032497 ± 0.0000000010) a.m.u. 

h = (1.054571628 ± 0.000000053) x lO'^^Js 

ke = (1.3806504 ±0.0000024) x IQ-^'^ JR-^ 

la.m.u. = (1.660538782 ± 0.000000083) x 10"^^^^ 

the value of the conversion factor for two ^He atoms should be taken as ttt- = (12.11928 ± 
0.00002) KA^. The major source of the uncertainty in this conversion factor is the Boltz- 
mann constant, followed by the Planck constant and the mass unit. 

The results of our calculations of the binding energy and effective range expansion pa- 
rameters for two He atoms are summarized in Table HI For conversion of energy between K 
and atomic units we use the value of the Boltzmann constant ks = 3.1668154 x 10^® a.u./ K 



We provide two types of results. The first group of calculations is performed with the ef- 
fective coupling constant used in previous three-body calculations, the second type is done 
with the best values recommended by NIST. The calculations are done in atomic units. 
For conversion between K and a.u. (the factor "c" discussed above) we also use the values 
recommended by NIST. This explains the difference between the values reported here and 
previously published results for the TTY potential. Evidently, the binding energy and the 

^ = 12.12 KA2 ^ = 12.11928 KA2 

m m 

Potential E2, niK a, A rg, A E2, inK a, A rg, A 

HFD-B(He) 1.6853 (0.4%) 88.60 (0.2%) 7.28 1.6921 88.43 7.28 

LM2M2 1.3034 (0.5%) 100.2 (0.2%) 7.33 1.3094 100.0 7.33 

TTY 1.3149 (0.5%) 99.82 (0.2%) 7.32 1.3210 99.59 7.33 

HFD-B3-FCII 1.5872 (0.4%) 91.19 (0.2%) 7.29 1.5938 91.00 7.29 

TABLE I. Binding energy and effective range parameters for different potential models calculated 
for exact and approximate effective coupling constants. Relative effect of rounding the coupling 
constant is given in parentheses. 

scattering length are quite sensitive to the often neglected small inaccuracy of the coupling 
constant. Even though the the coupling constant employed in many published calculations 
differs from the exact value only in the fifth figure, it is the third significant figure of the 
binding energy and the scattering length which is affected by this small inaccuracy. 

As the characteristics of near-threshold He-He states are quite sensitive to the accuracy 
of the effective coupling constant employed, it is useful to place physically reasonable limits 
on the accuracy of few-body calculations. Fig. [1] shows the coupling constant dependence of 
the He2 binding energy dependence for the TTY potential. The inset shows the variations 
of the two-body bound state energy within one-cr confidence interval for the effective cou- 
pling constant. This permits an estimation of the accuracy of the few-body near-threshold 
calculations at the level of 0.4 /xK, which corresponds to a relative error of the order of 
oc 3 X 10~^. This provides a natural limit to the accuracy of numerical calculations for this 
system, a limit deriving from the known precision of the fundamental constants. Explicitly, 
varying the effective coupling constant — from 12.11926 KA^ to 12.11930 KA^ we see the 
binding energy varying from —1.69231 mK to —1.69194 mK. Thus, numerical calculations 
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FIG. 1. He dimer binding energy as a function of the coupling constant. The inset shows the 
variation of the binding energy within the bounds determined by the inaccuracy of the fundamental 
constants. 

with an accuracy of more than four significant figures are not possible, since the potentials 
are defined in units of the Boltzmann constant ks- 



III. THREE-BODY OBSERVABLES 



The known strong correlations between the two-body scattering length and the energies of 
trimer bound states suggests that we can also expect sensitivity to the coupling constant in 
three-body calculations. The three-body results reported here were obtained using a general- 
purpose three-body code based on solving Faddeev equations numerically. The numerical 
technique used in this code is similar to the one reported in 22l, l25| . A very brief description 
of our approach to solving the Faddeev equations is given in the Appendix. The code is 
intended for public release and this work is a part of the release preparation. 



A. Three identical He atoms 

In order to determine the degree of sensitivity of three-body observables to small vari- 
ations of the coupling constant, we must first conduct a thorough analysis of the intrinsic 
numerical uncertainties and convergence properties when the coupling constant is held fixed 
at its recommended value. Convergence tables for bound states of the ^Hea trimer calculated 
with the LM2M2 potential are shown in Table HIl The bound state energies (with respect 
to the three-body break-up threshold) are presented for different numbers of grid points in 
spline solutions to the Faddeev equations. The number of grid points in the "cluster coordi- 
nate" X and in the "reaction coordinate" y are set equal, and the number of grid points in 
the angular coordinate z = (x, y) is varied independently. (The detailed description of the 
coordinate system can be found in [25]). The cut-off distances are set to R^; = 1200 a.u. for 
the cluster coordinate and Ry = 2000-1 /| a.u. for the reaction coordinate. Not surprisingly, 



the excited state is much less sensitive to the grid in angular coordinate: as is shown in 25| . 
the excited state is strongly dominated by the two-body s-state. 

To estimate the intrinsic numerical error we extrapolate the calculated values to an infinite 
grid and compare the calculated values with the extrapolated one. Estimating the error 
very conservatively we obtain the energy estimates E^ = (—4.0060 ± 0.0001) x 10~^ a.u.= 
-126.499 ± 0.003 mK for the ground state and E^ = (-7.21540 ± 0.0001) x 10^^ a.u.= 
-2.2784 ± 0.0003 mK for the excited state with the LM2M2 potential. Results obtained 
with other model potentials are given in Tab. IIIII 

We have also calculated the scattering length for a ^He atom scattered off the ■*He2 dimer. 



As our previous calculations have shown 22|, the box size is very critical for obtaining 



accurate results, and should be of order 3000 a.u. in the cluster coordinate. The results of 
scattering length calculations for the same set of potential models are shown in Tab. IIVI 

We can now address the issue of sensitivity to the coupling constant. How do small 
inaccuracies in the coupling constant used in many previous calculations affect the calculated 
observables? For this purpose we have introduced a coupling constant A = 12.11928/12.12 = 
0.9999406 in front of the potential and repeated the calculations using the most dense grid. 
In Tab |V] we compare results of such calculations with the results of previous calculations 
performed using similar techniques. Examining Tab |V] we find that the results obtained 
here confirm the results of earlier calculations performed with sparser meshes using a widely 



Ground state energy, a.u. 



N^ = Ny N^ = 1 N^ = 2 N^ = 4 N^ = 6 N^ = 8 oo 

20 -3.96752E-7 -3.99142E-7 -3.99137E-7 -3.99342E-7 -3.99468E-7 -3.99455e-07 

25 -3.98763E-7 -3.99602E-7 -4.00046E-7 -4.00415E-7 -4.00356E-7 -4.00360e-07 

30 -3.98560E-7 -3.99610E-7 -3.99967E-7 -4.00139E-7 -4.00111E-7 -4.00113e-07 

35 -3.98839E-7 -4.00237E-7 -4.00534E-7 -4.00670E-7 -4.00671E-7 -4.00671e-07 

40 -3.98781E-7 -4.00235E-7 -4.00476E-7 -4.00547E-7 -4.00562E-7 -4.00561e-07 

45 -3.98727E-7 -4.00230E-7 -4.00449E-7 -4.00541E-7 -4.00565E-7 -4.00563e-07 

50 -3.98764E-7 -4.00310E-7 -4.00513E-7 -4.00577E-7 -4.00607E-7 -4.00604e-07 

55 -3.98778E-7 -4.00273E-7 -4.00486E-7 -4.00571E-7 -4.00595E-7 -4.00593e-07 

60 -3.98753E-7 -4.00298E-7 -4.00494E-7 -4.00567E-7 -4.00597E-7 -4.00594e-07 

65 -3.98777E-7 -4.00288E-7 -4.00497E-7 -4.00577E-7 -4.00603E-7 -4.00601e-07 

70 -3.98763E-7 -4.00298E-7 -4.00496E-7 -4.00572E-7 -4.00602E-7 -4.00599e-07 
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Excited state energy, a.u. 


N, 


= Ny 


N, 


= 1 


Nz 


= 2 N, = 4: N, = 6 


N, 


= 8 


CX) 



20 -7.20149E-9 -7.21731E-9 -7.21671E-9 -7.21825E-9 -7.21908E-9 -7.21911e-09 

25 -7.20046E-9 -7.20577E-9 -7.20855E-9 -7.21090E-9 -7.21046E-9 -7.21047E-09 

30 -7.20196E-9 -7.20863E-9 -7.21089E-9 -7.21195E-9 -7.21172E-9 -7.21172E-09 

35 -7.20645E-9 -7.21541E-9 -7.21725E-9 -7.21809E-9 -7.21806E-9 -7.21806E-09 

40 -7.20136E-9 -7.21070E-9 -7.21216E-9 -7.21256E-9 -7.21262E-9 -7.21263E-09 

45 -7.20376E-9 -7.21344E-9 -7.21475E-9 -7.21528E-9 -7.21541E-9 -7.21541E-09 

50 -7.20411E-9 -7.21407E-9 -7.21527E-9 -7.21561E-9 -7.21578E-9 -7.21578E-09 

55 -7.20348E-9 -7.21311E-9 -7.21438E-9 -7.21487E-9 -7.21499E-9 -7.21500E-09 

60 -7.20384E-9 -7.21379E-9 -7.21495E-9 -7.21536E-9 -7.21555E-9 -7.21555E-09 

65 -7.20372E-9 -7.21345E-9 -7.21469E-9 -7.21515E-9 -7.21531E-9 -7.21531E-09 

70 -7.20374E-9 -7.21363E-9 -7.21481E-9 -7.21525E-9 -7.21543E-9 -7.21543E-09 

oo -7.20374E-9 -7.21358E-9 -7.21478E-9 -7.21522E-9 -7.21539E-9 -7.21540e-09 

TABLE II. Convergence tables for ^Hes bound state energies for LM2M2 potential. 



Potential Ground state energy Excited state energy 

£'3, a.u. Es, a.u. £'3, a.u. £3, a.u. 

Es/ks, niK E^/kB, niK 
(best grid) (extrapolated) [ (best grid) (extrapolated) 

HFD-B(He) -4.21430E-7 -4.21425E-7 -133.075 -8.68438E-9 -8.68438E-9 -2.74231 

LM2M2 -4.00602E-7 -4.00600E-7 -126.499 -7.21543E-9 -7.21540E-9 -2.27844 

TTY -4.00722E-7 -4.00720E-7 -126.537 -7.25704E-9 -7.25703E-9 -2.29159 



HFDBFCIl 


-4.15369E-7 


-4.15365E-7 


-131.163 -8.25861E-9 


-8.25861E-9 


-2.60786 


TABLE III. 


Bound state 


energies (with respect to the break-up threshold) 


for the 


^Hea trimer. 


Potential 


ai2, a.u. 


(best grid) 


Scattering length estimates 
ai2, a.u. (extrapolated) 




012, 


A 



HFD-B(He) 


230.416 


LM2M2 


218.051 


TTY 


218.936 


HFDBFCIl 


228.380 



230.424 
218.060 

218.945 
228.388 



121.93 
115.39 
115.86 
120.86 



TABLE IV. Scattering length for the ^He-^He2 atom-dimer scattering. 



adopted inexact value of the coupling constant. The values obtained for the TTY potential, 
however, differ from the previously reported values. As is mentioned above, in the previous 
calculations we used the value of the Boltzmann constant suggested by the authors of the 
potential. Here we are using the value recommended by NIST. It is also evident that the 
difference in three-body observables induced by the variation of the coupling constant at the 
level of 0.005% is well resolved by our code. We have summarized the results of calculations 
performed with the exact and the rounded value of the coupling constant in Table IVII 
It has been suggested by many authors that the system of three He atoms should demon- 



strate a nearly universal behavior 
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24l |. Using a separable interaction Pen'kov has 



shown 2M that all the low-energy parameters of the system of three ^He atoms can be 
described by a single dimensionless parameter within a few percent error. A similar obser- 
vation has been made by Braaten and Hammer [231] on the basis of effective field theory 
and by Platter and Phillips (27| who used an approach similar to 23], but used a higher 
order expansion and different regularization technique. In all these works the pair angu- 
lar momentum cut-off has been introduced and only the s-wave interaction has been taken 



Potential 



Scattering length estimates (rounded coupling constant) 



ai2, a.u. 



ai2, A 



012 



,A 



22] 



HFD-B(He) 


230.28 


LM2M2 


217.74 


TTY 


218.64 


HFDBFCIl 


228.21 



121.86 
115.22 
115.70 
120.76 



121.9 

115.4 

115.8 

n/a 



Potential 



Ground state energy (rounded coupling constant) 



£^3, a.u. 



Es/kB, niK 



E^/kB, mK 



M 



HFD-B(He) 
LM2M2 
TTY 
HFDBFCIl 



-4.21084E-7 
-4.00265E-7 
-4.00384E-7 
-4.15689E-7 



-132.968 
-126.394 
-126.431 
-131.264 



-132.98 
-126.41 
-126.40 
-131.26 



Potential 



Excited state energy (rounded coupling constant) 
El a.u. EykB, mK E^kB, mK [25 1 



HFD-B(He) 
LM2M2 
TTY 
HFDBFCIl 



-8.65867E-9 
-7.19200E-9 
-7.23353E-9 

-8.28580E-9 



-2.7342 
-2.2711 
-2.2842 
-2.6164 



-2.734 
-2.271 
-2.280 
-2.617 



TABLE V. Calculated observables with the coupling constant modified to match previous calcula- 
tions. 



into account. Therefore, when comparing results of direct full-configuration calculations 
with predictions based on universality arguments one can expect better agreement for the 
strongly s-wave dominated near-threshold states, with poorer agreement expected for the 
ground state of the trimer - which is not as dominated by the s-wave. Each of these theories 
has a parameter which should be fit to reproduce some three-body observable, and then 
predicts other three-body observables within a few percent error. A similar observation was 
also made by Delfino et al. 32|, who suggested that the ratio of neighboring bound state 
energies of an Efimov-like system is a universal function of the ratio of the dimer binding 
energy to the energy of the three-body ground state. We show this example in Fig. [2l The 
curve from the original work 32] has been recovered graphically, and the accuracy of this 
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^ = 12.12 KA^ ^ = 12.11928 KA^ 

Potential E^, niK £'3, mK a^, A E^, mK E^, mK ai2, A 

HFD-B(He) -132.968 (0.08%) -2.7342 (0.3%) 121.86 (0.06%) -133.075 -2.74231 121.93 

LM2M2 -126.394 (0.08%) -2.2711 (0.3%) 115.22 (0.15%) -126.499 -2.27844 115.39 

TTY -126.431 (0.08%) -2.2842 (0.3%) 115.70 (0.14%) -126.537 -2.29159 115.86 

HFD-B3-FCII -131.264 (0.08%) -2.6164 (0.3%) 120.76 (0.08%) -131.163 -2.60786 120.86 

TABLE VI. Energies of the bound states of '^Hes and the atom-dimer scattering length for different 
potential models calculated for exact and approximate effective coupling constants. In parethesis we 
show the relative difference between the results obtained with exact and inexact coupling constant. 

procedure is comparable with the observed discrepancy. It is, therefore, unclear, whether 



the small discrepancy comes from digitization or from the model employed in [32j. 

Similar correlations can be observed for the atom-dimer scattering length. As the atom- 
dimer scattering is dominated by the pole of the t-matrix corresponding to the near-threshold 
state of the trimer, the atom-dimer scattering length should behave as 

la 

VE2 - El [e* 7' 

where a is the two-body scattering length. Our numerical calculations confirm this simple 
observation extremely well. In Fig. |3]we show the ratio of the atom-dimer scattering length 
to the two-body scattering length as a function of the dimensionless parameter — ^= — . All 

the numerical results fall on a nearly perfect straight line. (Similar connection between the 
neutron-deuteron doublet scattering length and the energy of the three-body bound state is 



known in nuclear physics as the Phillips line 33|)- Again, we see that the shifts in scattering 



length induced by an inaccuracy of the coupling constant consistently follow the overall 
trend. 

Before turning to other examples of three-body systems that can be treated with our 
code, let us compare our results with other known full-configuration calculations. We have 
summarized the results reported for the ground state of the ^Hes trimer with the TTY 
potential in Tab. IVIII In the left column we summarize the results obtained in adiabatic 
hyperspherical (AH) representation, including single-channel 13|, l28| and multi-channel [ij, 
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calculations. In the right column the results of direct numerical solution of the 

11 
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FIG. 2. Ratio of the energies of two subsequent bound states of the Hes trimer as a function 
of the dinier binding energy. The dotted hne is recovered from Ref. [32] ( with digitization error 
comparable with the observed discrepancy). The difference between the results obtained with exact 
and inexact coupling constants is clearly resolved and consistent with the overall trend. 



Schrodinger [29| or Faddeev equations [18|, |26|, ISlI] are given. The result reported in J3l| was 
obtained with a restricted angular basis corresponding to the very simplest angular grid used 
in the present calculations, and agrees perfectly well with the value of -125.95 niK which 
we obtain in this simplified case. The overall agreement between the results of solving the 
equations directly is much better than those obtained within the AH approach. Although 
the AH approach provides an effective and reliable tool for studying few-body systems 
qualitatively, obtaining converged results is more difficult. The variations of the results due 
to inaccuracies in the conversion factors can not account for systematic underestimating 
of the binding energy which we see in adiabatic hyperspherical calculations. We want to 
emphasize that although obtaining converged results in the AH approach is technically 
more difficult, it should not be considered impossible. For instance, Suno and Esry [30] 
report E3 = —133.55 mK for the ground state and E^ = —2.7845 niK for the excited state 
energy of the trimer in the calculation with 35 AH channels using the SAPT2007 potential 
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FIG. 3. The atom-dimer scattering length (in units of the two-body scattering length) as a function 
of the distance to the pole in the three-body t-matrix. The difference between the results obtained 
with exact and inexact coupling constants is clearly resolved and consistent with the overall trend. 



AH Reference 



[JJ -106.1 niK (a) 
[28] -105 mK (a) 
[19] -125 mK (b) 
[14] -125 mK (b) 



Direct calculation reference 



[31] -125.9 mK 

[26] -126.39 mK 

[18J -126.40 mK 

[29] -126.2 mK 



-123.8±0.5 mK (b) Present -126.537 (126.431) mK (c) 

TABLE VII. Comparison of different results reported for ^Hes ground state with TTY potential, 
(a) Single-channel approximation has been used, (b) Full multi-channel calculation in AH repre- 
sentation, (c) The value corresponding to the rounded coupling constant is given in parenthesis 



10|. These values agree well with our result of E^ 



T33.589 mK for the ground and 



El = -2.78474 niK for the excited state. 
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IV. UNSYMMETRIC TRIMER 

Unlike the symmetric trimer, the trimer formed from two ^He atoms and one ^He atom 
has only one bound state. In order to approach the Efimov regime there should be two 
subsystems possessing large binary scattering lengths compared to the effective range. In 
the case of the unsymmetric trimer the ^He-^He scattering length is only twice as big as 
the effective range, and, therefore, the situation is quite far from the Efimov limit. The 
data published on the unsymmetric trimer is not as extensive as for the symmetric case. 
However, the ^He2^He trimer is a simple inhomogeneous three-body system, and as such 
it is an important example for the benchmarking of three-body codes. In this section we 
report tests similar to those shown above for the symmetric system. 

All the 2-body results are given in Table IVlnl As we can see from Table IVlnl the '^He-^He 



^ = 12.12 KA2 ^ = 12.11928 KA2 



Potential a, a.u. ro, a.u. a, a.u. tq, a.u. 

HFD-B(He) -34.382 18.40 -34.401 18.39 

LM2M2 -33.245 18.56 -33.263 18.56 

TTY -33.226 18.56 -33.244 18.55 

HFD-B3-FCII -34.050 18.45 -34.069 18.45 

TABLE VIII. Scattering length and the effective range for ^He-^He collisions. 

scattering length is less sensitive to the variations of the coupling constant, so the variations 
in three-body observables should mostly be due to the sensitivity of the ^He2 subsystem. The 
three-body data is summarized in Table HXl The sensitivity of the unsymmetric trimer bound 
state energy to small variations of the coupling constant is comparable to the sensitivity of 
the bosonic trimer excited state. The atom-dimer scattering length, however, is much less 
sensitive. 

Similarly to the case of the symmetric bosonic trimer, we have plotted the ratio of the 
^He-^He2 atom-dimer scattering length to the '^He-^He scattering length as a function of 
the dimensionless parameter —j^ — which characterizes the distance between the two-body 

threshold and the pole of the three-body t-matrix corresponding to the unsymmetric trimer 
bound state. All the data consistently follow the same trend (see Fig. H]). 
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^ = 12.12 KA^ 



^ = 12.11928 KA^ 



Potential E's, a.u. 012, a.u. i?3, a.u. 012, a.u. 

HFD-B(He) -5.3815E-8 (0.3%) 36.18 (0.03%) -5.3958E-8 36.17 

LM2M2 -4.5355E-8 (0.3%) 36.91 (0.03%) -4.5488E-8 36.90 

TTY -4.5270E-8 (0.3%) 37.23 (0.03%) -4.5404E-8 37.24 

HFD-B3-FCII -5.1442E-8 (0.3%) 36.57 (0.03%) -5.1582E-8 36.56 

TABLE IX. '^He2'^He bound state energy and ^He-^He2 atom-dimer scattering length. In paren- 
thesis we show the relative difference between the results obtained with exact and inexact couphng 
constant. 
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FIG. 4. The ^IIe-'^He2 atom-dimer scattering length (in units of the ^He-^He scattering length) 
as a function of the distance to the pole in the three-body t-matrix. The difference between the 
results obtained with exact and inexact coupling constants is clearly resolved and consistent with 
the overall trend. 

V. CONCLUSIONS 

We have performed highly accurate calculations of the system of three He atoms, including 
both symmetric and unsymmetric cases. We have studied the effects of small (0.006%) 
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variations of the interaction coupling constant on the energies of He dimer and trimer bound 
states and scattering lengths. Results of the calculations reported here are converged to 
5 significant figures for the most of the observables. All the results are consistent with 
previously published calculations and follow the same universal - potential independent - 
trend. We can summarize the sensitivity of the results to small inaccuracies in the coupling 
constant as follows: a 0.006% variation of the coupling constant induces 0.5% shift in the 
energy of the dimer, 0.3% shift in the energy of the excited state of the homogeneous trimer 
and the single bound state of the inhomogeneous trimer, 0.2% shift in the ^He-^He atom- 
atom scattering length, 0.06%-0.15% shift in the ^He-^He2 atom-dimer scattering length, 
0.08% shift in the energy of the ground state of the homogeneous trimer, and 0.03% shift 
in the scattering length for ^He-^He2 collisions. As the effective coupling constant for the 
potentials reported in units of temperature is known with the relative accuracy of 0.0002%, 
the physical limit for the accuracy of the binding energy of the He dimer, the excited state of 
homogeneous He trimer and the energy of the inhomogeneous He trimer is about 5 significant 
figures. This accuracy is achieved in the present calculations performed with a computer 
program being prepared for a public release. 

The authors hope that both the qualitative observations and the numerical results re- 
ported in this work can be used in the future for benchmarking various quantum few-body 
codes. 
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APPENDIX A 

Here we provide a very brief overview of our approach to solving the Faddeev equations 
numerically. More detailed and rigorous description of the equations and the numerical 
approach is being prepared as a separate publication. 



According to the Faddeev formalism [SJ] the wave function of three particles is expressed 
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in terms of Faddeev components $ 

^(xi,yi) = $l(Xi,yi) + $2(x2,y2) + $3(x3,y3) , 

where Xq, and y^ are Jacobi coordinates corresponding to the fixed pair a 



^mi^+m-y ' ^ P 1/1 



(2) 

Here r^ are the positions of the particles in the center-of-mass frame. The Faddeev compo- 
nents obey the set of three equations 

{Hq + 14(x„) - E)<l>„(x^, y„) = -V;(x„) Y. '^/3(x;3, y/3) 

/3^" , (3) 

a = 1,2,3 

where V^(x„) stands for the pairwise potential and Hq is the kinetic energy of the three 
particles. To make this system of equations suitable for numerical calculations one should 
take into account the symmetries of the physical system. As far as all the model potentials 
are central it is possible to factor out the degrees of freedom corresponding to the rotations 
of the whole cluster [35] . For the case of zero total angular momentum the reduced Faddeev 
equation reads 

(i^oa + Va{Xc) - E)^a{Xa, Vex, Z^) = -XaVaVaiXa) V] ^fi{Xl3, UlS, Zjj) . (4) 

Here i/o« = -fe - || - (^ + i)^(l - .^)^^, 

Xa |Xa| ; 

Va = Wal , (5) 

_ (Xg, Ya) 
Za ) 

XaVa 

and the coordinate transformations between different system of Jacobi coordinates follow 
from the definition of the Jacobi coordinates (El) . 



The asymptotic boundary condition for bound states consists of two terms [34] 



^{x,y,z)r^ (t)2{x)e ^^y + A{-,z 



X 



-fc3(a;^+J/2)^ 



1 } 

14 



where 02 (a^) is the two-body bound state wave function, ky = \/E2 — E^, k^ = y/—E^, E2 is 
the energy of the two-body bound state and E-i, is the energy of the three-body system. The 
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second term - which corresponds to virtual decay of the three body bound state into three 
free particles - decreases much faster than the first one which corresponds to virtual decay 
into a particle and a two-body cluster. In our calculations we neglect the second term in 
the asymptotic introducing the following approximate boundary conditions for the Faddeev 
component at sufficiently large distances R^ and Ry 



d^Hx,y,z)[y=R^ ^ ^ ^^) 

^{x,y,z)[y=R^ ^' 

In order to solve the equations numerically we introduce a basis of Hermit splines sat- 
isfying the boundary conditions ([6]) and use orthogonal collocations to calculate a discrete 
matrix analog of the Faddeev operator 



More detailed description of the discretization procedure can be found in [25j (see also [36 1 
where we describe the procedure of constructing optimal non- uniform grids automatically). 
We solve the system of linear equations iteratively exploiting factorability of the left-hand 
side of Eqs. @ for preconditioning. In particular, we introduce localized components 

Ta = {Hoa + Vj^ - S^E)^^ (7) 

which - due to the asymtotic properties of the Faddeev components - has much better 
spacial localization than the original Faddeev component. The equations for the localized 
component read 

Ta = —VaiCapiHojB + VflS*^ — Sj^E)^ T/J + Ca-yiHo^ + V^S^ — S^E)^ T^) (8) 

The transformation ([7]) makes it possible to reduce the rank of the linear system essentially, 
^ically by a factor of 3. Application of Eqs. (jH]) to atom-dimer scattering can be found in 



22|. A similar idea has been discussed in [37|. 
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